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Abstract 

An amplifier has been designed for optimal use of position sensitive 
thermal neutron detectors using the principle of resistive charge division. 
The important points in this optimization are: high counting rates and 
good spatial resolution. This amplifier is built as a hybrid circuit and 
is now used on several new instruments at the ILL. It consists of a fast 
low noise current pre-amplifier, a gaussian shaping circuit based on a 
4th order active filter and an essentially noiseless baseline reconstruction. 
In this paper, we present a rather complete theoretical analysis of the 
problem that lead us to the choices made above, and allows for an optimal 
adaptation to other situations. An analysis of unwanted, secondary effects 
is also worked out. 



1 Introduction. 

The original motivation for the design of this electronics was given by a project 
called SANS-2MHZ, in the framework of the Millennium program at the Institut 
Laue Langevin, the large instrument renewal effort the ILL is making in order 
to be competitive in the 21st century. The aim of the project was to make 
a large area neutron detector with a counting rate of at least 2MHz at 10% 
dead time correction. One of the possibilities considered was to have an array 
of independent, 1 meter long, 1-dim position sensitive proportional detectors. 
These detectors are thin (8 mm diameter) and are based upon resistive charge 
division. We will present a general analysis of the problem of spatial resolution 
as a function of the electronic noise sources and take this as a guiding principle 
to explain our design. 

We will also analyze unwanted secondary effects. Four different contributions 
are considered: extra passive resistance, the presence of a blocking capacitor, the 
finite input impedance of the amplifiers and the presence of a certain capacitive 
load on both sides. 
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2 Noise analysis of the resistive charge division 
problem. 

Position sensitive detection based upon charge division is a very old and well- 
known technique, described in, for example, PP, or in [2]. In a classical publica- 
tion by Radeka [3], it is shown that the spatial resolution that can be obtained 
by resistive charge division along an RC line is determined by the temperature 
T, the capacity C and the collected charge Q s in the detector: 



Ws 

when using perfect amplifiers and a signal shaping that is optimal with respect 
to this problem. However, the fundamental hypothesis formulated to obtain 
that result is that the generated charge along the electrode is fast compared 
to the time scale of RC. For thermal neutron detection in cylindrical position 
sensitive detectors, using a wire and He-3 gas conversion, this hypothesis is 
not satisfied at all. Indeed, a typical detector has, say, 6KH,/m resistance and 
lQpF/m capacity. The total charge collection time is not negligible compared to 
the time constant, which is about 60 ns ; on the contrary, these 60 ns are almost 
insignificant compared to the charge collection time. This charge collection time 
is of the order of hundreds of nanoseconds, and is essentially determined by the 
difference in time of arrival between the first and the last primary charges along 
the tritium and proton tracks arriving in the amplification region on one side, 
and the finite drift time of the avalanche ions in the high field region, which 
generates the current signal, on the other side. 

Anticipating an (evident) result which will follow in our analysis, the spatial 
resolution will improve with a higher signal (collected total charge) level ; so we 
will have to work with shaping functions which are broad (in the time domain) 
compared to the time constant of the RC line in order to collect most of the 
useful signal. This, on one hand, invalidates the result in equation^m our case, 
but on the other hand, simplifies strongly the analysis ; indeed, we do not have 
to consider the dynamic behavior of an RC line, with all its Bessel responses 
and so on 0] . In the low- frequency limit where we are have to work in order to 
collect our charges, the detector essentially behaves as a lumped resistor. 

The preamplifier noise contribution can always be modelled as a series volt- 
age noise source and a parallel current noise source. The voltage v(t) and the 
current i(t) of both sources are random processes which are supposed to be 
stationary, and gaussian. They are hence described by their spectral noise den- 
sities. They could be correlated because they are obtained as equivalent sources 
of the noise sources internal to the preamplifier, but this correlation will be 
neglected. When the first stage of the preamplifier consists of an operational 
amplifier, the constructor of the circuit gives those spectra in the data sheets. 
For all practical purposes, they are a constant (it is white noise). Only in the 
case of a MOSFET entry, the 1/f noise could have a significant contribution. 
We will suppose in a first approach that the amplifier is perfect if we put the 
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Figure 1 : Equivalent diagram describing the noise sources in the charge division 
problem in the ideal case. 



noise sources outside. 

The final element which will complete our noise analysis is a feedback resis- 
tor going from the output to the input. The Johnson noise of that resistor will 
add quadratically to the parallel noise source of the amplifier. We have now 
completed the description of our basic noise circuit in resistive charge division 
which is shown in figure ^ Although the amplifiers are modelled as operational 
amplifiers (and in our case, we really do use such amplifiers), our analysis ap- 
plies to just any feedback amplifier structure. All sources are considered to be 
independent of each other. The measured quantities are the currents in Rfbi 
and Rfb2- We consider the amplifiers ampl and amp2 to be perfect. Node 1 
is hence kept at V\ and node 2 at vi. Solving for the current in the feedback 
resistors, we have: 



h = Hhi +h — Hh + 



h = Hh 2 +T-2 + ith + 



1>2 - Vl 

"anode 
^-anode 



(2) 



Let us consider that i\ includes ith t and that ii includes ith 2 - When calcu- 
lating the sum and the difference of both signals: 



/- = h — Ii = ii — 12 
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(3) 
(4) 



1+ = h+h=h + *2 
The spectral densities of the signals combine in square and the cross terms drop 
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out, so we find, using the Johnson noise expression for the resistor noise sources: 

J 2 - 2i 2 ■ + 16fcT 4- o ^noise (t -\ 
L - - /z noisc + d , + 8 j?2 W 
'"-anode ""-anode 

4 = 2 moise ( 6 ) 

Here, it is understood that * no j se = ^amp + ^fcT/i?^, iamp being the equivalent 
current noise density of the amplifier, and « no i se is the equivalent voltage noise 
density of the amplifier. Although these expressions arc in general spectral 
quantities, in the adopted hypothesis of white noise, they reduce to numbers. 
In order for the thermal noise to be the dominant contribution, R ano de should 
be in the range: 

v 2 ■ RhT 
noise „ t) „ ohjL in\ 

-%jjfr < Ranode < (<> 

noise 

If we take as an example an amplifier and feedback current noise of 2.0pA/y/Hz 
and a voltage noise of InV/ V Hz, this interval goes from 125S1 to 8KQ. 

If the time-invariant transfer function from input current to output voltage 
is H(f), then the root mean square (r.m.s.) of the output noise is given by: 

vlm.s. = / il\H{f)\ 2 df (8) 
Jf=o 

In our case of white noise, this reduces to: 




Vr.m.s. =in\l ^ jH(f)\ 2 df (9) 

So knowing the input current noise spectral density, we have to multiply by a 
number which depends only on the overall transfer function (shaping function) 
to find the 1 sigma deviation on the output voltage. Of course, this operation 
also applies to the sum and the difference signals. The contribution of the 
electronic noise on the spatial resolution of the detector can be determined as 
follows. Consider the dimcnsionless parameter p = (A — B)/(A + B). Here, A 
and B are instantaneous samples of the output signals of the amplifiers at their 
peaking value. In the ideal case, p runs from -1 to 1 when the physical position 
runs from to length L. The error on p due to the noise, considering that the 
noise on A-B and A+B is not correlated (as can easily be verified assuming 
identical transfer functions), is given by: 



+ v 2 



/ 2 2 

SrmsP = 1 (10) 



This translates in a FWHM position resolution of: 



v ■ +r.m.s. + u -r.m.s. 

5 FWRM X = v2in2L (n) 

^signal 
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Here, v s jg na i is the sum of both voltage signals for an "average" event. 

Let us now analyze the general behavior of the influence of the shaping 
function on the spatial resolution. Imagine an initial charge q, which is deposited 
instantaneously at the entrance of the amplifier, so i(t) = qS(t). The output of 
the amplifier is nothing else but the impulse response of the system times q, and 
H(f) is the Fourier transform of that impulse response. Of course, an overall 
amplification factor doesn't influence the spatial resolution as the signal and the 
noise are multiplied by it, so in order to compare transfer functions, one should 
normalize them such that the peak value (the maximum of the impulse response) 
is always equal to 1 . It is easy to work out that the amplitude normalized impulse 
response scales in time as: 

hr(t) = h^t/r) (12) 

and: 

H T (f)=rH 1 (rf) (13) 

In these expressions, the subscript 1 corresponds to a time-normalized impulse 
response, and r corresponds to the time scale parameter (the shaping time 
constant). 

The input noise calculations are independent of the transfer function, and 
are related to the output noise by the factor given in equation [5] For an am- 
plitude normalized transfer function, this factor then gives the absolute "noise 
performance" of a shaping function, so we will call this factor the absolute noise 
performance factor (ANPF) of the shaping function. It should be as low as 
possible. It is easy to determine that the ANPF scales as yfr when using the 
expression ! 131 So we have the at first sight remarkable conclusion that the faster 
the amplifier (the shorter the shaping time) the better the position resolution. 

If we can define a time scale parameter which is independent of the shape of 
the impulse response, then we can define a 'relative noise performance factor' 
by dividing its ANPF by y/r. We then have a quantity which can compare 
different shapes of transfer functions according to their merit. The problem 
with this definition is of course that a judicious choice of time scale parameter 
can artificially give an advantage to one or another shape so one should be 
prudent in drawing conclusions. 

We can define the gain G of the amplifier as the maximum of the impulse 
response. Writing H(f) the normalized transfer function, the full transfer func- 
tion of an amplifier is then GH(f). In that case, f s ig na i in formula |lll is equal to 
GQ S . As expected, this factor G drops out of the expression for the resolution. 



W2i 2 . (l+p 2 )+ B 16fcT + p , nmse 

, _ V n01SC anode Vode 

%WHM* = V2hi2L^RNPF-Y (14) 

Expression^Jexpresses the resolution due to electronic noise as a function of the 
shape and time scale of the overall transfer function, the event charge produced 
by the detector, the anode resistance and the equivalent voltage and current 
noise sources of the preamplifier. The hypothesis is put forward of an otherwise 
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ideal amplifier and full dynamics. If the sensitive region occupies only a fraction 
D (0 < D < 1) of the whole dynamics of the quantity p = (A - B)/(A + B), 
then the resolution will have to be scaled up with an extra factor equal to 1/D, 
see section [5] As long as D is not too far from 1, we can neglect the influence 
of these extra components in the equivalent circuit on the noise calculation. 



3 Gaussian Shaping Function. 

Although the pure pole approximation to a gaussian filter is a well-known item, 
we shortly sketch its derivation in order to make clear all our normalization 
conventions. The ideal gaussian shape as an impulse response would be propor- 
tional to: 

exp(-i 2 /2) 

/2T 



^perfect W = Jkz ( 15 > 



its associated transfer function (in the variable cu = 2irf ) equals exp(— cu 2 /2), 
or, in the Laplace variable s, exp(s 2 /2). Of course, as such, it is not realizable. 
Using the well-known approximation scheme of this transfer function into a 
stable, pure-pole filter by expanding D(s) = l/(H(s)H(—s)) as a polynomial up 
to order 2n, and constructing H(s) = l/£)+(s), where D+(s) is the polynomial 
of order n with the same zeros as D(s) in the left half complex s-plane, we find 
for n = 4: 

4.90... 

V-D+Cs) = 4.90... + 11.42.. 7s + 10.86.. .s 2 + 5.07.. .s 3 + s 4 ^ 

where we took the value of the numerator such that H(0) = 1 (because in the 
process, the amplitude normalization became arbitrary). The impulse response 
(the inverse Laplace transform of H(s)) approaches the perfect gaussian response 
(up to a time shift) and corresponds to a realizable, stable pure pole filter of 
4th order. We arbitrarily define its time scale r to be equal to 2ir, which 
corresponds to about the time where the response is "significantly different from 
0" (eg, the total visible width of the pulse). 1/r also corresponds closely to the 
cut-off frequency of this low pass filter. The final thing to do is to normalize 
the amplitude of the impulse response to 1. The peaking time of the impulse 
response occurs at 0.3114r, and it turns out that we have to divide the above 
transfer function by 2.507 in order to obtain this amplitude normalization. So 
we obtain finally: 

h T (t) = 3.86 e - 8 - 52t / T cos^-3.86e- 7 - 42t / T cos^ (17) 

r r 

+ 35.33e- 8 - 52 ^ sin ^ - 10.29e- 7 - 42 ^ sin ^ 

T T 

The corresponding transfer function is: 

3046. It 

T ^ S ' ~ 7635.2 + 2833.4sr + 429.0s 2 r 2 + 31.9s 3 r 3 + s 4 r 4 (1 ' 
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Now that we have the normalized transfer function of this filter, we can calculate 
the ANPF (by using H(s — i2nf) ) and it turns out to be: 

ANPF = 0.3788^ (19) 

Filling this into the expression for the spatial resolution, we obtain: 



J 2j Lisc( 1 +P 2 )+«^ + -? 2i ^ 
w noise "anode anode 

*FWHM* = 0.446LV7^ (20) 

Let us take as an example a 1KVL resistance, and the same noise sources as 
cited earlier. If we have a detector that delivers a harmonic average of OApC 
per event and we apply a gaussian shaping time of lAfis, then on a 1 m long 
detector we should be able to achieve a resolution of 12 mm along the detector. 
If the resistance is 6.5KQ, the resolution in the middle of the detector is 5.6 
mm, and on the borders 6.7 mm. The numbers in this last example have not 
been chosen at random: they correspond to an experimental situation we have 
set up. Measurements using a narrow neutron beam on a prototype Reuter 
Stokes detector with an overall resistance of 6.5Ktt lead to an observed FWHM 
resolution 5.7 mm in the middle of the detector and a deterioration towards 
7 mm near the ends of the detector, confirming our calculation. In this mea- 
surement, wide word lengths (10 bits) were used to code the position. If we 
measure the r.m.s. noise of the sum and difference output signals of gaussian 
amplifiers having a gain of 3V/pC in this setup, we measure 6.0mV rm s for the 
difference and 4.1mV rms for the sum signal. Calculations give 5.7mV rm s for the 
difference signal, and 3.8mV rms for the sum signal, showing an agreement with 
the measurements with less than 10% deviation. 

Let us compare this gaussian transfer function with the "standard" RC-CR 
function: its amplitude normalized impulse response is given by : 

KccAt) = -e 1 "- (21) 

r 

where r is its peaking time and its transfer function equals: 

H rccr (s) - ^^—f (22) 

The ANPF of this RC-CR filter equals e/(2V2)Vr = 0.961^. However, in 
order to really compare the ANPF values, we should adapt the r in this defini- 
tion to the one we used in the gaussian filter. We will propose identical peaking 
times: in that case, we have to replace t by t — 0.311441t' in the above ANPF 
expression, leading to ANPF = O^^. So we observe that the gaussian 
filter, for identical peaking times, improves the spatial resolution by about J^2% 
over the result one would obtain using the standard RC-CR filter. 
If we take a rectangular shaping function with width r, 

/isquare(t) = 1 - step(t - r) (23) 
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the transfer function is: 



-^square (s) 



l-e~ 



TS 



(24) 



-s 




and the ANPF equals \/t/2 — 0.707y^r. Again, in order to compare the r 

values, we take the width of the square function to be the FWHM, the peaking 
time not being applicable in this case, of the gaussian approximation which 
equals 0.38171r. So we should fill in r = 0.38171r' in the expression of the 
ANPF, to obtain that the ANPF = 0.4369^- Again, the gaussian filter 
improves the spatial resolution by 15% over the result one would obtain with a 
square filter for identical FWHM times. 

All this indicates that the gaussian transfer function as derived above is a 
very acceptable solution to the shaping problem for resistive charge division. 

4 The electronic implementation. 

The analogue part of the circuit implemented at the I.L.L. in the frame of the 
SANS-2MHZ program takes on the following form: a pre-amplifier which acts 
as a current amplifier, a gaussian filter implemented as an active filter and a 
baseline correction circuit. 

4.1 The preamplifier 

The preamplifier needs to satisfy 3 criteria: the amplifier noise contributions 
have to be such that they don't deteriorate significantly the spatial resolution, 
the input impedance has to be as low as possible and it has to stand very 
high counting rates. The fact that there is anyway a resistive noise source (the 
resistive anode) eliminates the noise advantages a charge amplifier has. Indeed, 
in order to stand very high counting rates, we opted for a current amplifier, 
which suffers no pile up effects. The transfer function of the circuit is chosen to 
be a first order system with a time constant of the order of the time constant 
of the detector, about 60 ns. In order to achieve a very low input impedance, 
a high open loop gain is required, as explained in subsection 15.31 The high 
bandwidth combined with the high open loop gain made us decide to go for a 
bipolar operational amplifier. 

Expression 1201 is the guiding principle to select the operational amplifier, 
which has to have a very high open loop unity gain bandwidth. The best values 
one can obtain for the current noise are around 1.5pA/y/Hz, which doesn't 
contribute significantly to the resolution as long as the anode resistance is below 
~ lOKtt. This comforts our choice of a bipolar circuit, such as the CLC425 from 
National Semiconductor, unfortunately now not in production anymore. The 
feedback resistor doesn't have a significant influence compared to the current 
noise of the amplifier as long as it is above about 12Kfl. A feedback capacitor of 
about 5 pF lets us obtain the right bandwidth in this case, which is a satisfying 
value. 
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Figure 2: Basic scheme of the preamplifier. 



The basic scheme that we obtain for the structure of the preamplifier is 
shown in figure El 

As long as the voltage noise is a few nV/V~H~z, and the anode resistance is 
not too low (above ~ lKfl) this term doesn't contribute significantly to the 
spatial resolution. If a resistive protection circuit against accidental discharges 
is required, care should be taken to include it in the feedback loop (as shown in 
figure [2J, otherwise its resistance adds to the input impedance. One should be 
warned, however, that the resistive part of such a circuit will add to the voltage 
noise of the circuit. 

In order to obtain a larger gain, a T-type feedback structure can be used 
(see figure EJ. An extra advantage of using a current amplifier over a charge 
amplifier is that no pole-zero cancellation network is needed. The preamplifier 
is AC coupled to the gaussian filter. 

4.2 The gaussian filter. 

The fact that for neutron detection, time constants of the order of the microsec- 
ond are used allows us to use operational amplifiers (designed for the telecom 
industry) in active filters to implement the gaussian transfer function. Fac- 
torizing expression 1181 into two second-order contributions, we notice that the 
Q-factors of both circuits are quite low and that hence a simple Multiple Feed- 
back Network |S] can be used to implement them. Although possible, passive 
LC filters require rather large self inductance values on this time scale, so we 
reserve that kind of filters to faster circuits. 
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Figure 3: The baseline correction circuit: extraction of the baseline. 



4.3 The baseline correction circuit. 

Given the fact that the gaussian impulse response is unipolar, there is a problem 
of a shifting baseline at high counting rates. A simple clamping circuit has 
the disadvantage of clamping on the noise, so we would have a noisy baseline 
correction. In order to solve this problem, we use the circuit displayed in figure 
[3] In that figure, i? m translates the input voltage into a current. If the voltage 
is positive, the feedback loop is closed through resistor R2 and diode D2 ; the 
voltage at the output is 0. If the input voltage is negative, the loop is closed 
through Rl and Dl. The output voltage is now — Rl / R- m V^ n , which is a positive 
quantity. Taking i? m = Rl, if the signal from the gaussian shaper is applied to 
this circuit, then the output equals minus the baseline which is negative. It is 
now sufficient to low pass filter this output and add it to the input signal with a 
simple active adder circuit to restore the baseline. The path R2/D2 prevents the 
feedback loop from ever opening completely, and hence keeps the operational 
amplifier from saturating. If the application is continuous irradiation, then 
the baseline will vary very smoothly and a simple first order low pass (RC) 
circuit with a cutoff frequency below say \KHz can do. However, in pulsed 
time of flight applications, the baseline also fluctuates more rapidly, and a more 
sophisticated filter might be necessary in order to filter all noise and glitches, 
but reacting fast enough to the changing baseline. 

4.4 The conversion circuitry. 

An analogue peak detector with threshold on the sum signal triggers 2 single shot 
12 bit ADC converters with a dynamic range of 5 V. The 11 most significant 
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Figure 4: Equivalent circuit of all secondary effects we are considering. 



bits of each conversion are used to compose a 22 bit word which serves as 
an address in a position lookup table. This lookup table is materialized by 
an EPROM circuit ; the result of the lookup operation is the position of the 
impact in an 8-bit word. The quantization noise of the ADC (= 5V/2 11 l/>/T2 = 
0.7mVrms ) should be smaller than the noise at the output of the amplifiers. 
It should be noted that here, the right choice of the gain to map the signal 
range to the dynamic range, is important. The quantization noise added to 
the resolution by truncating the result to an 8 bit word results in a FWHM 
resolution of L/2 8 l/vT2\/81n2 = 0.00266L, which should be added m square 
to the resolution obtained in 1201 Taking our measurements again, which had a 
spatial resolution of 5.7 mm using wide word lengths, when coding on 8 bits, 
using the above formula, we should increase theoretically the resolution to 6.3 
mm. Measurements lead to 6.6 mm of resolution when coding on 8 bits. 

5 Secondary effects. 

In order to perform the above analysis, we used a "perfect" circuit. In reality, 
there are unavoidable secondary effects which have an influence on the perfor- 
mance. We will study the 4 most prominent ones, and make the assumption 
that the situation remains symmetrical. First we will study the presence of an 
additional passive resistance. Interesting in its own right, it also introduces the 
notion of "equivalent resistance" for any other contribution to the reduction of 
the dynamics, and hence makes it possible to compare the importance of several 
different contributions. A most important influence is played by the blocking 
capacitor, which is a necessary component to protect the electronics from the 
high bias potential applied to the resistive electrode. It has two effects: one is 
a reduction of dynamics (just as the extra resistance), the other is to introduce 
an undershoot in the signal. The finite input impedance of the preamplifier 
can also contribute significantly to the loss in dynamics. Finally, the presence 
of a capacitive load (due to the wiring up of the detector) will not influence 
the dynamics, but will introduce a noise contribution that can deteriorate the 
spatial resolution. 

We will study all these effects one by one but we will not study their combined 
effects. The overall scheme is displayed in figure 0] The ideal case would give us 
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a relationship between the position x and the collected charges qi and q 2 given 
as follows: 

Qi = (l- i)<lm><l2 = 2<lm ( 25 ) 
We can hence introduce an ideal dimensionless position p as follows: 

<?2 - qi n - r 2 

P = = (26) 

qi +q2 n+ r 2 

When p = — 1, we hit position x = 0, and when p = +1, we hit position x = L. 
We will observe that several effects change this relationship into: 

p=^^ = D r -l^l (27) 
qi +q2 n+ r 2 

where D is a number between and 1, and is called the dynamical factor. 
5.1 Additional resistance. 

Imagine that the resistive electrode continues to be resistive beyond the sen- 
sitive area, or that some resistance is present in the conducting circuit before 
reaching the — supposed perfect — current amplifier at impedance (virtual 
ground). Let us suppose that the case is symmetrical and that an extra resis- 
tance ^p ara site 1S added at both sides. We have now, for an injection at point 
x of a total length L: 

* = ^ ; 2 + + ; r parasite ™ 

11 zl parasite 

v = ?in ;\ + ;r asite ^ 

- rL " r Al parasite 

If we now calculate the range that can be taken up by p = (q 2 — qi)/{qi + q 2 )- 

p= n-r 2 = R r 1 ~r 2 

R + 2 r parasite ^ + 2r parasite ri + r2 

we observe that the dynamical factor is: 

D= - ~ ^/parasite 

R + 2rp arag ^ e R 

So the calibration, (q 2 — qi)/(qi + q 2 ) as a function of (ri — r 2 )/(ri + r 2 ), 
stays linear, but the slope is reduced. The second approximation is valid if the 
parasitic resistance is small compared to the useful resistance of the resistive 
electrode. 
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5.2 Blocking capacitor. 

If we introduce a blocking capacitor C^\ QC , or shortly C, on each side (see figure 
U) we find as current transfer functions (injected current to collected current): 

• f \ 1 + Cr 2 s 

11{S) = 2 + C(r 1+ r 2 )s ^ 
■ I \ l + Cr lS 

12{S) = 2 + C(r 1+ r 2 )s (33) 
If we look at the time domain responses of these transfer functions: 

2t 

■ U\ e c(-i+-2) (n ~ r 2 ) r 2 

H{t) = C(r 1+ r 2 )i + WTV 2 m (M) 

m - e ^f^ ( V l) + -^hm (35) 

C(ri+r 2 y n+r 2 

The time signals clearly contain two parts: the second term is exactly what we 
would have if we were in the ideal case, and the first term is the effect of the 
capacitors. It is an exponential decay with a time constant of C(r\ + r 2 )/2. 
It has moreover the opposite amplitude behavior than the ideal term: indeed, 
when the injection point goes more to the left, so that the i\ current increases, 
the extra term becomes more and more negative. On the other hand, when the 
injection point goes to the right, so that the ideal term almost vanishes, the 
extra term becomes positive and more and more important. 

The two extra terms in i\ and i 2 are opposite, so that they cancel in the sum. 
This is something important: on the sum signal, the effect of the extra capacitors 
vanishes. This means that the sum signal contains the timing information, 
independent of the position of the impact. 

We should of course convolve this injected current with the overall impulse 
response of the shaping amplifier. There are two contributions: 

h*h(t) = (n-r 2 )V{t) + - h(t) (36) 

i 2 *h(t) = (r 2 - ri )V(t) + — 5— (37) 

n + r 2 

In these expressions, we introduced the function V(t): 

V(t) = ^-*h(t) (38) 

V(t) doesn't depend on the injection point along the electrode, but only on 
R = ri + r 2 . We notice two things: while the second term in each expression 
gives exactly what we wanted to obtain, the first one disturbs the response 
in a non-proportional way. This means that the individual maxima of the two 
functions will in general not occur together in time, and the relative deformation 
of the signal will depend on the injection point. If we would take those individual 
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maxima as charges to be collected, the calibration curve would be non-linear. 
However, if we can take a common sample time T for both curves which is 
independent of the injection point, and which corresponds to the maximum of 
the impulse response (G), we find the following: 



Qi = {n-r 2 )V{T) + ■= — G (39) 

n + r 2 

q 2 = (r 2 -r 1 )V(T) + —^-G (40) 



= g 2 -gi = n - T2 L _ 2(n +r 2 )V(T) 
g2 + qi n +r 2 \ G 

So we see here that the calibration curve corresponds to a straight line, with a 
reduction factor of the dynamics that is equal to: 

D = 1 _^ = 1 _ 2 _i_ (e -^ (()) ^ (42) 

When we compare the expression of the loss in dynamics with the one for a re- 
sistor at both ends, we can identify both formulas in the approximation of small 
deviation from unity. In that case, it is as if the coupling capacitor introduced 
an equivalent parasitic resistance: 

req=^(e- 2 ^*M*)) t=T <«> 

If we assume that the time constant CR is larger than the time scale (and of 
course the sampling time) of the shaping function, then the exponential with 
which we have to convolve looks essentially like a unit step function. If we 
replace it in the convolution integral, we have: 

req -dc l- h{t)dt (44) 

The equivalent resistance of the coupling capacitor turns out to be simply in- 
versely proportional to the capacity times a constant that only depends on the 
shaping function. It really has a meaning of an equivalent resistance. It turns 
out that for a Gaussian 4th order filter function, this expression is numerically 
equal to: 

r cq ~ 0.1576^ (45) 

Here, r is the time scale of the Gaussian shaper as we defined it. 

A remark of a technical nature is maybe due: the capacitance is of course the 
dynamic capacitance (small signals) for the capacitor under the static bias load 
that one has in mind. For a perfect capacitor, that doesn't make any difference, 
of course ; however, for ceramic capacitors using the X7R dielectric, for example, 
that dynamic capacity is reduced by 40% when working at the nominal voltage 
I 
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A second problem, introduced by the presence of blocking capacitors, is the 
fact that there can be an undershoot on the "large" signal. This undershoot 
is a problem in the case a baseline correction is used. We can estimate an 
upper boundary of the value of the undershoot as follows: given that we are 
just outside of the time lapse where the shaper function is important (if it is 
a well-peaked function such as a Gaussian), the undershoot is essentially the 
sample value of the "non-ideal" part of the impulse response at that moment. 
Of course, there will be a slight contribution of the shaper response, so the 
actual undershoot will be slightly less than this value. But we will obtain a 
reliable upper boundary. Let us calculate this value now for n = and r 2 = R: 

ui ~ RV(T') (46) 

We sample at t = T', assuming that h(t) has no contribution anymore there. 
This can be worked out as follows: 



1 f T ' 
h(t) = -^e-T^T' / h{t)e^dt (47) 
™^ Jt=o 



nT' 

e " rc ' 



The last integral can be extended to infinity (because h is supposed to be neg- 
ligible outside of the range 0-T'). As H(s) normally has poles with a real value 
much more negative than the exponential coefficient, this operation is allowed. 
We can replace the integral by a sample point of the Laplace transform of h 
(the shaper transfer function). In that case, we finally obtain as our undershoot 
estimation (in absolute value): 

^"W'-* rE {-w) <48) 

To obtain the relative undershoot, we have to divide that expression by the gain 
of the shaper function. Knowing the scaling relations of the transfer function 
in the case of unity gain, we see that the undershoot only depends on the 
ratio of T" /(RC). In the case of a gaussian filter, then taking T' = t gives us 
numerically: 

3046. lxe- 2 * 

Ul ~ 7635.3 + 2833.4(2z) + 429.01(2;r) 2 + 31.9(2a;) 3 + (2a;) 4 ~ ( } 

where: 

x = « 1 (50) 

The undershoot introduces an essentially random error in the case of a baseline 
correction, with an amplitude limited to the above amount. It depends on the 
image that is projected onto the detector. It should hence be smaller than the 
relative error introduced by the noise over the full dynamics of the signals. 

We use 68nF capacitors as blocking capacitors in our application. With 
1.4/xs gaussian shaping, this corresponds to an equivalent resistance of 3.217 
which is really negligible. The maximal undershoot, given a dynamics of 5V, 
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and a detector resistance of 6.5K£l, is 6.3mV. This is of the same order of 
magnitude as the noise levels (see our previously calculated example). We see 
here that it is in fact the undershoot which gives us the most severe condition 
on the value of the blocking capacitor. 



5.3 Input impedance. 

The input impedance of the preamplifier circuit also has as an effect to reduce 
the effective dynamics of the charge division. We will give an estimation of this 
effect. The input impedance of a current amplifier with a transimpedance Z(s) 
determined by a passive feedback loop, is easily found out to be equal to: 

Z m (s) = (51) 

where A(s) is the open loop transfer function of the amplifier. We will take it 
that Z(s) is represented by a first order system. In figure this is the feedback 
resistor times the gain factor T of the T-network, and the time constant 
ipre is given by the feedback resistor and its parallel capacitor. It is a bit more 
delicate to write down the open loop transfer function A(s) of the operational 
amplifier: indeed, large bandwidth amplifiers usually have very complicated and 
not very well known transfer functions. However, because of the gaussian filter, 
we only use a limited part of that bandwidth, and we will assume that we can 
model the open loop gain by a first order system with dc gain A i and time 
constant t ;. The impedance then becomes: 



T fr R ib( 1 + t ol s ) 
A ol (l + t we s) 



Z m(°) = , 7 V (52) 



If we now define the following equivalent elements: 

^fb^fb 



r = 



(53) 



= t oF (54) 



we can rewrite the input impedance as: 

r + Is 
1 + tpres 



Z in(*) = i ,T ( 55 ) 



At first sight, one would be inclined to neglect the pole in this expression. 
Indeed, the preamplifier is supposed to represent a current amplifier, and hence 
be much faster than the gaussian shaping filter that follows. In that case, the 
input impedance of the amplifier is equivalent to the series circuit of a resistor of 
value r and an inductor of value I. We tried this, and to our surprise the results 
are numerically quite different from the results obtained when we don't neglect 
this pole. So this pole cannot, in the cases we're interested in, be neglected. 
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If we analyze the structure of the input impedance, we can easily show jS] 
that this corresponds to the passive circuit displayed in figure[5] Using the same 
approach as in the case of the blocking capacitor, we find that the currents take 
on the form: 

( ) _ (n - r 2 )(l - rt pro ) - ^^re * , 1 + ^pre g(m) 
1[) (2? + (r 1+ r 2 )Ve) 2 + 21 + (n + r 2 )t pre 6{W 

_ (n - r 2 )(l - rtp re ) - ^T^? t 1 + n t we ( ] , } 

Again, the sum of both signals is not affected by the presence of this input 
impedance. We next convolve these expressions with the transfer impedance of 
the shaping amplifier, which in our case is the 4th order gaussian filter function. 
Finally, sampling at the time when the sum signal attains its maximum, we 
obtain the sampled charges qi and q 2 at each side. The calibration remains 
linear, but a reduction of dynamics is introduced by this input impedance. 
This allows us to define an equivalent resistance. The result of these algebraic 
operations is: 

rRt wc V{a)+l{2r + R-RV{a)) 
r «l ul U Rt W c + 2rt W e(l-V(a)) + 2lV(a) [ ™> 

In this expression, we defined R = r\ + r 2 the anode resistance, and: 

2r + R 

a = 2Z + ^ pr e TgaUSS (59) 
The function V(a) is now defined as follows: 

/•0.3114 

V(a) = ae'- ' 3114 "' / h^e^dt (60) 

using the time normalized gaussian filter function hi(t). We observe that the 
equivalent resistance doesn't depend on the injection point (is no function of r\ 
or r 2 individually), which confirms the linearity of the calibration, but it is a 
function of the total anode resistance R. Working out numerically the integral, 
we find: 

_ a(3046.1e-°- 3114a + (a - 12.02)(153.3 - 19.9a + a 2 )) 
( °°~ (76.8 - 17.0a + a 2 )(99.5 - 14.8a + a 2 ) ( ' 

We can obtain an interesting result when we take the limit of big R, meaning, 
much bigger than 2Z/tp r e or r. a then converges to a limiting value given by: 

Tgauss , . 

Q lim = -7 (62) 

tpre 

and the expression for the equivalent resistance reduces to: 

lim K 1 - ^( a lim)) , t r / \ , aQ s 
r equi = ^ prc + rV >hm) (63) 
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We now observe that for large enough detector resistance R, the equivalent 
input resistance of the circuit is a well-defined quantity. In cases where the 
input impedance is a crucial parameter, we should optimize the ratio of the 
gaussian shaper and the preamplifier time constant such that the associated 
equivalent resistance is minimized. In fact, the perfect match occurs when the 
preamplifier time constant is equal to l/r = t Q \. Usually, this time constant is 
too big to be used (it should still be smaller than the gaussian shaping time), so 
a less optimal choice has to be made. Diminishing the gain of the preamplifier 
also reduces the input impedance. 

As an example, we have taken a preamplifier based on a CLC425 opera- 
tional amplifier (modelled with open loop dc gain of 20000 and corner frequency 
lOOfc-ffz), a feedback resistor of 12K £1 in parallel with a 5.6pF feedback capaci- 
tor, and a T-factor of 13.0. With a 1.4/is gaussian shaping, the above calculation 
leads to an equivalent resistance of 26.90 when the detector is 500ft ; we have 
measured 26f2. The limiting value is 21.0O. 



5.4 Capacitive load. 

We will now analyze the effect of an extra load capacitor on both sides of 
the detector, as displayed in figure 01 This corresponds partly to the detector 
capacity itself, but mainly to the extra capacitive load that is introduced by the 
conductors leading from the sensitive electrode up to the input of the amplifiers. 
An order of magnitude is 1 pF per cm of (coaxial) connection wire. If we add 
this capacity C on both sides of the noise equivalent network in figure ^ the 
noise currents now become: 

h = Hhi + H ~ Hh + „ sCvi 

-"-anode 

h = ith 2 +i2+ Hh + -=r — - sCv 2 (64) 

anode 

The spectral densities of the sum and difference signals then take on the follow- 
ing form: 

I'M) = Koise + ^T + Koise^f^ + J^) ^ 
anode anode 

'+(/) = Koi S e + Koise- 2 f 2c2 ( 66 ) 
So to both densities, we add the quantity: 

&(/) = 8< oise7 r 2 / 2 C 2 (67) 

This time it is clear that these noise spectral densities are a function of frequency 
f. We cannot apply the simple factor anymore as in equation E| In fact, our 
noise formulas are of the form: a + bf 2 , so we can still use for the a-part, but 
we now need a second number for the b-part. The extra part equals: 

/>oo 

vl m . s .-c = ^ 2 v 2 nQise C 2 / f\H{f)\ 2 df (68) 

Jf=0 
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We define now the capacitive noise number rj: 



V = / f\H(f)\ 2 df (69) 

It can easily be verified that r\ scales as 1/r, and not, as is the case for the white 
noise contribution, as r. So this capacitive load noise becomes more and more 
important when the time scale becomes small. We can call ^/rj the capacitive 
noise performance factor (CNPF) if we work with the amplitude normalized 
transfer function. For the gaussian 4th order approximation, we can work out 
that: 

CNPF.^ (70) 
V T 

This means that we have an effective extra contribution to the sum and differ- 
ence noise current densities (considered as white noise), equal to: 

ff _ 2V27rCNPF %Qise C 
lc ~ ANPF [71 > 

For the gaussian shaping function, this turns out to be: 



■eff . 



6-70 %gigeg (72) 



As an example, a l.OnV/vHz voltage noise density over a capacity of 100 pF 
(1 meter of coaxial cable) and a time constant of 1/j.s gives us an effective 
noise current of O.&lpA/ \J H z. This is still smaller than the other noise current 
contributions (but its effect should start to become visible). 

In order to add this effect to the spatial resolution, it is sufficient to replace 
the sum and difference currents by their combination, squared, with this effective 
current. 



6 Conclusion. 

The principles leading to the design of the front end electronics for resistive 
charge division in thermal neutron detection in high counting rate applications 
are exposed. The choice we made for the analogue part consists of three sub 
functions: a preamplifier which acts as a current amplifier, based upon a fast, 
bipolar operational amplifier ; a 4th order gaussian shaping circuit based upon 
simple active filter circuits, and a noiseless baseline correction. A detailed anal- 
ysis is presented leading to the expression feauation ll4|) which gives the spatial 
resolution as a function of the detector signal, the noise sources of the pream- 
plifier and the gaussian shaping time constant. A comparison of the gaussian 
transfer function with the standard RC-CR and the square pulse serves as a 
motivation for the gaussian shape. A short description of the practical imple- 
mentation as hybrid circuits follows. 
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A profound analysis of several secondary effects are given. Two important 
effects are due to the blocking capacitor. Its equivalent dynamical resistance is 
given in equation 1451 and an upper boundary to the undershoot it introduces 
is given in equation Also a detailed analysis of the influence of the input 
impedance of the amplifier is presented, leading to indications how to optimize 
this quantity if necessary. Finally, the noise contribution due to the capacitive 
load is worked out in equation [721 

Some experimental results confirm the calculations presented here. 

Although our analysis was focussed on the design of the front end electronics 
for the SANS-2MHz project at the Institut Laue Langevin, the results are of 
a sufficiently general nature to allow them to be a guide for other designs in 
position sensitive detection. 
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